A NEW CEMENT TO GLUE NON-CONFORMING GRIDS WITH 
ROBIN INTERFACE CONDITIONS: THE FINITE ELEMENT CASE 
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Abstract. We design and analyze a new non-conforming domain decomposition method based 
on Schwarz type approaches that allows for the use of Robin interface conditions on non-conforming 
grids. The method is proven to be well posed, and the iterative solver to converge. The error analysis 
is performed in 2D piecewise polynomials of low and high order and extended in 3D for Pi elements. 
Numerical results in 2D illustrate the new method. 
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1. Introduction. Our goal in writing this paper is to propose and analyze a 
non-conforming domain decomposition generalization to P.L. Lions initial idea |30| 
in view of an extension of the approach to optimized interface conditions algorithms. 
This type of algorithm has proven indeed to be an efficient approach to domain de- 
composition methods in the case of conforming approximations |14j , [24j . This paper 
presents the basic material related to so called zero**^ order (non optimized) method 
in case of finite element discretizations, see [18] for a short presentation. In the com- 
panion paper [2] , the case of the finite volume discretization has been introduced and 
analyzed. 

We first consider the problem at the continuous level: Find u such that 

= / in f7 (1.1) 
C{u) ^ g ondn (1.2) 

where C and C are partial differential equations. The original Schwarz algorithm is 
based on a decomposition of the domain into overlapping subdomains and the reso- 
lution of Dirichlct boundary value problems in the subdomains. It has been proposed 
in [30j to use more general boundary conditions for the problems on the subdomains 
in order to use a non-overlapping decomposition of the domain. The convergence rate 
is also dramatically increased. 

More precisely, let fl he a C^'^ (or convex polygon in 2D or polyhedron in 3D) 
domain of M'^, d = 2 or 3; we assume it is decomposed into K non-overlapping 
subdomains: 

n = uf=in'. (1.3) 

We suppose that the subdomains fl'^, 1 < k < K are either C^'^ or convex polygons in 
2D or polyhedrons in 3D. We assume also that this decomposition is geometrically con- 
forming in the sense that the intersection of the closure of two different subdomains, 
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if not empty, is either a common vertex, a common edge, or a common face in 3D. Let 
rifc be the outward normal from fl'^ . Let {Bk/)i<k,e<K,k^i be the chosen transmission 
conditions on the interface between the subdomains (e.g. Bk,e — + What we 
shaU caU here a Schwarz type method for the problem p.ip - (|1.2p is its reformulation: 
Find {uk)i<k<K such that 

C{uu) - / in 17'= 
C(ufc) = g on d^t n dVt 

BkAuk) = BkAui) on dn^ n dn^. 

Let us focus first on the interface conditions Bkj- The convergence rate of as- 
sociated Schwarz-type domain decomposition methods is very sensitive to the choice 
of these transmission conditions. The use of exact artificial (also called absorbing) 
boundary conditions as interface conditions leads to an optimal number of iterations, 
see [33], [32] , [H] and [H]. Indeed, for a domain decomposed into K strips, the num- 
ber of iterations is if, see [32j . Let us remark that this result is rather surprising since 
exact absorbing conditions refer usually to truncation of infinite domains rather than 
interface conditions in domain decomposition. Nevertheless, this approach has some 
drawbacks: 

1. the explicit form of these boundary conditions is known only for constant 
coefficient operators and simple geometries, 

2. these boundary conditions are pseudo-differential. The cost per iteration 
is high since the corresponding discretization matrix is not sparse for the 
unknowns on the boundaries of the subdomains. 

For this reason, it is usually preferred to use partial differential approximations to the 
exact absorbing boundary conditions. This approximation problem is classical in the 
field of computation on unbounded domains since the seminal paper of Engquist and 
Majda [17j . The approximations correspond to "low frequency" approximations of 
the exact absorbing boundary conditions. In domain decomposition methods, many 
authors have used them for wave propagation problems [15], [16], [29], [6], [13], [27] 
and ^ and in fluid dynamics [3T], [5D]. Instead of using "low frequency" in space 
approximations to the exact absorbing boundary conditions, it has been proposed to 
design approximations which minimize the convergence rate of the algorithm. Such 
optimization of the transmission conditions for the performance of the algorithm was 
done in, [24] , [25| for a convection-diffusion equation, where coefficients in second order 
transmission conditions where optimized. These approximations are quite different 
from the "low frequency" approximations and increase dramatically the convergence 
rate of the method. 

When the grids are conforming, the implementation of such interface conditions on 
the discretized problem is not too difficult. On the other hand, using non-conforming 
grids is very appealing since their use allows for parallel generation of meshes, for local 
adaptive meshes and fast and independent solvers. The mortar element method, first 
introduced in [5], enables the use of non-conforming grids. It is also well suited to 
the use of the so-called " Dirichlet-Neumann" ([2^) or "Neumann-Neumann" precon- 
ditioned conjugate gradient method applied to the Schur complement matrix ([28], 
[3], |35j). In the context of finite volume discretizations, it was proposed in [34J to use 
a mortar type method with arbitrary interface conditions. To our knowledge, such 
an approach has not been extended to a finite element discretization. Moreover, the 
approach we present here is different and simpler. 
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Our final project is to use interface conditions such as 002 interface conditions 
(see [23], [IS]), this will be developed in a future paper. Here we consider only interface 
conditions of order : Bk.e = + Q^fc- The approach we propose and study was 
introduced in [18] and independently implemented in [26j for the Maxwell equations 
but without numerical analysis. 

2. Definition of the method and of the iterative solver.. We consider the 
following problem : Find u such that 

{Id-A)u = f mfl (2.1) 
u^O on dfl, (2.2) 

where / is given in L^(n). 

The variational statement of the problem (|2.ip - (|2.2p consists in writing the problem 
as follows : Find u £ 11^(^1) such that 

{WuWv + uv)dx ^ I fvdx, VveH^in). (2.3) 



Making use of the domain decomposition p.3p . the problem (|2.3p can be written as 
follows : Find u E i?o(f^) such that 



K ^ K 



Let us introduce the space i?^(il'^) defined by 

Hlin'') = {ipe H\n''), = over dnndn''}. 

It is standard to note that the space i?o(ri) can then be identified with the subspace 
of the if-tuple v= (ui, vk) that are continuous on the interfaces: 

K 

V = {v^{vi,...,VK)el[Hlin''), 

k=l 

\fk,e,k^e, l<k,e< K, Vk = ve over dfl'' n dn'^}. 
This leads to introduce also the notation of the interfaces of two adjacent subdomains 

j.kJ ^ g^k ^ g^e 

In what follows, for the sake of simplicity, the only fact to refer to a pair {k,£) 
preassumes that F*^'^ is not empty. The problem (|2.3p is then equivalent to the 
following one : Find ugV such that 

/ (VufeVwfc + UkVk) dx = fkVkdx, Vv G V. 

k^i-l^" k=i-^^'' 

Lemma 1. For v £ Yl^^i Hl(Q''), the constraint Vk — ve across the interface T^'^ 
is equivalent to 

K 

\fp EE {pk) e Yl H-^''^{d^^) wtth pk = -pe over F'^-^ 
fc=i 



K 



_H--i/2(ao'=) < Pfc:'f^fc >//i/2(ao'=) - 0- (2.4) 



k=l 
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Proof The proof is similar to the one of proposition III. 1.1 in 12J but can't be 
directly derived from this proposition. Let p = {pk) G IlfcLi H^^^^{dil'') with pk = 

—pi over r'^'^, in {H^l^^ {T^'^))' sense. Then, there exists over each Vl^ a lifting of 
the normal trace pk in i?(div, fi'^). The global function P, which restriction to each 
^l^ is defined as being equal to the lifting, belongs to 7?(div, fi) and is such that 
(P.n)|gf^fc — Pk- Let now v_£V. From the previously quoted identification, we know 
that there exists v S H}^{^) such that V|s^fc — Vk- In addition, 

vdwP - / PVv = 0. 
n Jo. 

On the other hand, 

/ -vdivV - [ PVv = V( / vdwP - / PVv) 

so that (|2.4p is satisfied. 

Reciprocally, let v = {vi,...,vk) & Ylk=i^*i^'^) ^^^^ that (|2.4p is satisfied. Let 

X G r*^'^, and let 7a; C 7^; C Fj; C f j; C F*^'^ be open sets. There exists a function ip in 

1 /2 

V{Tx) such that (p{y) — 1 for all y in ■fx- With any q e (-ffoo (-l^a;))'' associate 
p = (pk) defined by 

H-^/^dn") <Pk,Wk >HU^[Qn'') = (ifi/^(r,))' < 9'<^'^fe Hl'^'^ir^)^'^'^'' ^ H^^^id^l''), 

and pj = 0, Vj 7^ fc,^. 

Then, by construction, p £ HfeLi H^^^^{dil'^) and = — over T'^'^ . Hence from 
(EH), 

K 

H-i/2(afifc) <Pk,Vk >//i/2(9Qfc)= 0. 

fe=i 

We derive 

_f/-i/2(an'=) <Pk,Vk >ffi/2(anfc)= -H-i/2(aaf) <pe,ve >H^/^{dn')^ 

thus, 

(<^(r.))' < >ffor(r.)=(<^r.))'< ><^(r.)' 

1/2 

and this is true for any q G (^00 i^x))' , hence Lpvk — fvi over F^,, and thus 

Vk = ve over j^, G F*''^. 
We derive Vk = vi a.e. over F*^'^, which ends the proof of lemma [TJ 
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The constrained space is then defined as follows 

/ K \ / K 



vu = vi and qt = ~qi over T^'^} 
and problem (|2.3|) is equivalent to the following one : Find (m, p) e V such that 



k=i'^^'' k=i 

fkVkdx. 



K 

E 

k=l 



Being equivalent with the original problem, where pk — over dVl'^ (remind that / 

is assumed to be in L'^{Vl) so that actually belongs to H^^^'^{dVt^)), this problem 
is naturally well posed. This can also be directly derived from the proof of an inf-sup 
condition that follows from the arguments developed hereafter for the analysis of the 
iterative procedure. First, let us describe this algorithm in the continuous case, and 
then in the non conforming discrete case. In both cases, we prove the convergence of 
the algorithm towards the solution of the problem. 

2.1. Continuous case. Let us consider the interface conditions of order : 

Pk + auk = -pt + aui over F'"'^ 
where a is a given positive real number. 

We introduce the following algorithm : let {uk,Pk) G Hl{^^) x be 
an approximation of (u,p) in fj'^' at step n. Then, (u)!^^, p)!^^) is the solution in 
HliQ}') X iJ-i/2(5f7fc) of 



{Vul'^^^Vk + ul+^vk) dx- H-i/2(QQ.>') <Pl'^^,Vk >Hi/2(an'=) 



= / fkVkdx, yvkeHlin'^) (2.5) 
Jni' 

< pI+^ + aul+\vk >rK^=< -p? + aul\ >^-.^ Vz;,. G i/oo'(F'='0 (2-6) 

It is obvious to remark that this series of equations results in uncoupled problems set 
on every . Recalling that / e i^(ri), the strong formulation is indeed that 

- Au^+i + ul+^ = h over 9.^ 
^ ' ^ul+^ = -ri' + au'l over F*^-^ 



= over (2.7) 

OXlk 



From this strong formulation it is straightforward to derive by induction that if each 
p^, fc = 1, ...,K, is chosen in J^^ iJ^/^(F*'''^), then, for each k, I < k < K , and n > 
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the solution belongs to H^{fl'') and belongs to JJ^, H^/'^{V^'^) by standard 
trace results (p)!^^ — — p" + a(it" — mJ!^^)). This regularity assumption on p\ will be 
done hereafter. 

We can prove now that the algorithm (|2.5p - (|2.6p converges for all / G L^(ri). As 
the equations are linear, we can take / = 0. We prove the convergence in the sense 
that, in this case, the associated sequence ((u)J,p)J))„ satisfies 

lim (ll"fellifi(fi'=) + bfcllH-i/2(an'=)) =0, for 1 < fc < if. 

We proceed as in { [30j . [T^ ) by using an energy estimate that we derive by taking 
v^. = u^^^^ in (12. 5p and the use of the regularity property that pj!^^ G L'^{dVt^) 

that can also be written 

JQ'' 

f (t n+1 I n+l\2 / n+1 n+l\2\ j 

e 

By using the interface conditions (|2.6p we obtain 

= ^E/ + (2.8) 

Let us now introduce two new quantities defined at each step n : 

K 



and 

K 



fe=i e^k • 

By summing up the estimates ()2.8p over fc = 1, iiT, we have, 
so that, by summing up these inequalities, now over n, we obtain 



n=l 



We thus have lim„ >oo -B" = 0. Relation (|2.7p then implies : 

1™ bfcllH-i/2(9ofc) = 0, for k^l, K. 
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which ends the proof of the convergence of the continuous algorithm. 

Theorem 1. Assume that f is in L'^{Q) and {p^l)i<k<K G Yl^H^/'^iT^'^) . Then, 
the algorithm I12.5\) - I[27S\) converges in the sense that 

J^oo ^""^ ~ ^^'^^H^n") + \\Pk - PkWH-^/^dn")) ^0, forl<k<K, 

where Uk is the restriction to fl^ of the solution u to k2.1\) - W7^) . and pk ~ over 
dn^, 1 <k< K. 

2.2. Discrete case. We introduce now the discrete spaces. Each il'^ is provided 
with its own mesh 7^*^, 1 < k < K, such that 

For T S T/^, let hr be the diameter of T {hr — sup^ y^rp d{x, y)) and h the discretiza- 
tion parameter 

h = max hk, with hk — max hx- 

l<k<K TeT^ 

At the price of (even) more techniques and care in the forthcomming analysis, possible 
large variations in the norms of the solution u\Q_k can be compensated by tuning of 
hk- This requires that the uniform h is not used but all the analysis is performed 
with /ife. For the sake of readability we prefer to use h instead of /ifc. Let pT be the 
diameter of the circle (in 2D) or sphere (in 3D) inscribed in T, then ctt = ^ is a 

measure of the nondegeneracy of T . We suppose that Tf^ is uniformly regular: there 
exists a and t independent of h such that 

VT e , (TT < cr and rh < hr- 

We consider that the sets belonging to the meshes are of simplicial type (triangles or 
tetrahedron), but the analysis made hereafter can be applied as well for quadrangular 
or hexahedral meshes. Let ViviiT) denote the space of all polynomials defined over T 
of total degree less than or equal to M. The finite elements are of lagrangian type, of 
class C". Then, we define over each subdomain two conforming spaces Yj^ and by 

Yl: - {vh,k e C"(n'), Vh,k\T e Vm{T), VT e T,^}, 

^ {vhM e Yj^, f''i,fc|an'=nao = 

In what follows we assume that the mesh is designed by taking into account the 
geometry of the F*^'^ in the sense that, the space of traces over each F*^'^ of elements 
of YI^ is a finite element space denoted by 3^^'^. Again, at the price of more notations, 
this assumption can be relaxed. Let k be given, the space is then the product 
space of the 3^^'' over each £ such that F*^'^ ^ 0. With each such interface we associate 
a subspace W^'^ of 3^,^'^ in the same spirit as in the mortar element method [8] in 
2D or [5] and TT] in 3D. To be more specific, let us recall the situation in 2D. If the 
space consist of continuous piecewise polynomials of degree < M, then it is readily 
noticed that the restriction of to F'^^^ consists in finite element functions adapted 
to the (possibly curved) side F'^'^ of piecewise polynomials of degree < M . This side 
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has two end points that we denote as x^'^ and that belong to the set of vertices 
of the corresponding triangulation of V^'^ : Xq'^, x\"^ , x'^'^^^x'^'^ . The space Wj^'^ is 
then the subspace of those elements of y^'^ that are polynomials of degree < M — 1 
over both [a;o' ,Xi' ] and [x'^'_i,x'^'^]. As before, the space Wj^ is the product space 
of the W^'^ over each £ such that r'^'^ ^ 0. 
The discrete constrained space is then defined as 



{{Ph,k + auh.jc) - {-Ph.e + oiUhi))iph,k,i = 0, yiph,k,e e Wj^'^}, (2.9) 



and the discrete problem is the following one : Find [u^tP^ G Vh such that 

'yvh = {vh,i,-vh,K) € nf=i^^ 

/ i'^Uh.k'Vvh.k + Uh,kVh.k) dx / Ph,kVh.kds = / fkVh.kdx. (2.10) 

We introduce the discrete algorithm : let {u]^k>Phk) ^ -^h ^ ^ discrete ap- 

proximation of {u,p) in fi'^ at step n. Then, {Uh\^>Ph\^) ^he solution in Xj^ x 
of 

(v?/;^+iV«^,fe + K\W,k) dx- [ Pl\\h.kds = / fkVh,kdx, \fvh,k G (2.11) 



(PhJc^ + aKjc^)'^h,kd = I {-Pl,e + aulj,)^/jh,k,e, 'y'^h,k,i & W^;". (2.12) 



In order to analyze the convergence of this iterative scheme, we have to precise the 
norms that can be used on the Lagrange multipliers p^. For any p g Hs^i L^{di^''), 
in addition to the natural norm, we can define two better suited norms as follows 



llPll-i, 



K K 

EE 

k=i e=i 
i^k 



V 



\Pk\ 



where 11. 11 _i stands for the dual norm of iJnn(r'^'^) and 



K 



\k=l 



3 (90*=) 



We also need a stability result for the Lagrange multipliers, and refer to [3] in 2D 
and to the annex in 3D, in which it is proven that. 

Lemma 2. There exists a constant c» such that, for any Ph,k.e ^h'^ > there 
exists an element w^'^'^ in that vanishes over dVl^ \ T^^^ and satisfies 



Ph,k4W^-^''^ > \\ph.kA\^ -1 



(2.13) 



with a bounded norm 

..h,k,e 



\\w^''''\\HHn^)<c4Ph,kA^-i,^,,- (2.14) 



Let T^k^i denote the orthogonal projection operator from onto W^'^ . 

Then, for v e Trk,e{v) is the unique element of VF^ ' such that 

(nkAv) - v)^ = 0, e Wf/. (2.15) 

We are now in a position to prove the convergence of the iterative scheme 
Theorem 2. Let us assume that ah < c, for some constant c small enough. 

Then, the discrete problem Ii2.10\} has a unique solution (uinP^) G Vh- The algorithm 

i2.11\) - P2.1}3fj is well posed and converges in the sense that 

liiji^ I \\K,k-Uh,k\\m{n>')+'^\\pl,k,e~Ph,kA\ -i =0, for 1 < k < K. 

i^k * ^ ' 7 



n >oo 



Proof. For the sake of convenience, we drop out the index h in what follows. We 
first assume that problems (|2.10p and (|2.1ip - (|2.12p are well posed and proceed as in 
the continuous case and assume that f ~ 0. From (|2.15p we have 

and (|2.12p also reads 

pI+^ + a7rfe,,«+i) = iTkA~P7 + aw") over F'^'^. 
By taking Vk = u)!^^ in (|2.1ip . we thus have 

/ (iv<+ip+iurip)dx 

Then, by using the interface conditions (|2.12p we obtain 

= J- E / i'^k,£{p7 - auDfds. 

It is straightforward to note that 

(TTkAPi - au'Dfds < f {p'l - au'lfds 
(p" - m:tk{u^) + a7r£^fc(u") - au^fds 



{p'l - a7r^,fc(u^))^ + a^(7r£,fc(u^) - u'lyds 
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since {Id — tt^ is orthogonal to any element in Wf^^ . We then recall that (see 

[5] in 2D and [5] or [U equation (5.1) in 3D) 

^(^,,fc«) - u'l fds < c/i||u,"||^v2(r^..) (2.16) 

< c/l||u"||^i(j^(!^. 

With similar notations as those introduced in the continuous case, we deduce 

and we conclude as in the continuous case: if cah < 1 then lim„^oo E"" = 0. The 
convergence of towards in the norm follows. The convergence of p^! in the 
^-5^pfc/^ norm is then derived from ((27T3)l and ((2lT|) . Note that by having / = 
and (u",p") = prove that (it"+^,p"+^) — from which we derive that the square 
problem (|2.11|) - (|2.12|) is uniquely solvable hence well posed. Similarly, having / = 
and getting rid of the superscripts n and n + 1 in the previous proof gives (with 
obvious notations) : 

E + B < cahE + B. 
The well posedness of (|2.10p then results with similar arguments. 

We shall address, in what follows this question through a more direct proof. This 
will allow in turn, to provide some analysis of the approximation properties of this 
scheme. 

3. Numerical Analysis.. 

3.1. Well posedness.. The first step in this error analysis is to prove the sta- 
bility of the discrete problem and thus its well posedness. Let us introduce over 
(nf=i Hl{n^) X nf^i L^idfl'')) X nf=i HUn'^) the bilinear form 

A' „ ^ r 

a((u,p), w)) = / (VufcVwfc -|-MfeWfe)(ia; - / PkVkds. (3.1) 
The space Hl{VL^) is endowed with the norm 

Lemma 3. 

There exists a constant /3 > independent of h such that 

K 

k=l 

H{uh,Pj,Vh)) > miAhW* + ll£jUi,JkJU- (3.2) 
Moreover, we have the continuity argument : there exists a constant c > such that 

K 

k=l 

«((2ih,Pj,«J) < c(||^,||* + ||£jUi)(||l!,|U). (3.3) 
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Proof of lemma [3} In (16) and (17), we have introduced local HQ{Tk i) functions 
that can be put together in order to provide an element Wh of IlfcLi -^h ^^^^ satisfies 



K 

V / pkWkds > \\p P 



Let us now choose a real number 7, < 7 < ^ (where c, is introduced in (|2.14p ) and 
choose Vh — ILh ~ IMLh (EH) that yields 



a{.{uhiP,X^h)) = (VufeV(ufe - jwk) + Uk{uk - jwk)) dx 

K 

- / Pk{uk~ lWk)ds 
k^Jd^" 

As already noticed in section [2?2l we can write 

PkUkds = / {{Pk + omk,i{uk)Y - {pk - aTrk^e{uk))'^)ds 
1 

4a 
1 

4a 
1 

4a 



T- / ((7i'fe,£(-P£ + aui,)Y - {pk - aiTk,t{uk)f )ds 
" /r'=.« 

- T~ / ((P^ ~ ""^)^ ^ (^''t ^ a7^fe,Hwfc))^)rfs 

- 7~ I {{Pe - ct^e,k{ue) f ~ {pk - aTTk,e{uk)f)ds 



1 
4a 



a {-Ki,k{ut) - ue) ds 



so that 



K ^ K 

Uk - '^k,eiuk)yds < cah\\ui-^\\ 



/ PkUkds <jY ( 

Going back to (|3.4p yields 

d{{Uh,P,):llh) > (1 - cah)\\uf^\\l - l\\uh\\*\\wh\\* + tIIP/JI^ 1 
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>(i-ca;i)|K||2+7||pJp_,^^_i:||^J|2 
>(i-ca;i)||u,||2 + (7-2M)||pj|2_,^^. 



(3.4) 



Due to the choice of 7, we know that, for ah small enough, (j3.2p holds. The continuity 
(|3.3p follows from standard arguments (note that the norm on the right hand side of 
(|3.3p is not the ||.||_i ^-norm), which ends the proof of lemma [3l 



From this lemma, we have the following result : 

Theorem 3. Under the hypothesis of theorem\^ the discrete problem k2.10fl has 
a unique solution {ui^,p^) S Vh, and there exists a constant c > such that 

IhhW* + llP/JI-i,* ^ c||/||l2(o). 
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From lemma 131 we are also in position to state that for any (u^,j3^) G V/i, 

\\u-Uh\\* + ll£-£/JI-i,* ^ c(||u- w,J|» + ||p-£^ILi) (3.5) 
and we are naturally led to the analysis of the best fit of {u,p) by elements in Vh- 

3.2. Analysis of the best fit in 2D. In this part we analyze the best approx- 
imation of {u,p) by elements in Vh- As the proof is very technical for the analysis of 
the best fit, we restrict ourselves in this section to the complete analysis of the 2D 
situation. The extension to 3D first order approximation is postponed to the next 
subsection. 

The first step in the analysis is to prove the following lemma 
Lemma 4. Assume the degree of the finite element approximation M < 13, there 
exists two constants Ci > and C2 > independent of h such that for all ri{_k in 
yi,k p H}^{Y^'^), there exists an element ipg^k in W^''^ , such that 

im.k + T^k,d'ni,k))lpi,k > Cl\\rii^k\\l2(-pk,e), (3.6) 

ll'0€,fc||L2(r'=.f) < C2|l?7<!,fe||L2(r'=.o- (3-7) 

Then, we can prove the following interpolation estimates : 

Theorem 4. For any u G H'^{n) n H^{fl), such that Uk = ■ti|o'c G H'^+'"'{n''), 
1 < k < K with M ~ 1 > m > 0, u— {uk)i<k<K o-nd let pkj = over each T^'^ - 
Then there exists an element Uy^ in Y\^^iX^ and p^ = (pkih), Pkth G W^'^ such that 
{iLhTPh) saiis/j/ the coupling condition \2.9jl . and 

k=l k<e 
WPkih -PkAlfj-ii^^k.i) ^ cQ;ft.^+™(||ufc||H2+™(o'--) + \\ue\\H^+<^{ni)) 

+ ch'^n\Pk,\\^i.^^^.,) 

where c is a constant independent of h and a. If we assume more regularity on the 
normal derivatives on the interfaces, we have 

Theorem 5. Let u G H'^{n) n Hl{n), such that Uk = u\^k G H'^+'^ {D.'') , 

1 < k < K with M — 1 > m> 0, u= {uk)i<k<K and pk^t — is in Hi'^"^(rk/)- 

Then there exists U/^ in Y[^=i-^h '^"■'^ P/j ~ {Pkih), Pkth G W^^'^ such that {Uf^,p^J 
satisfy i2.9\) . and 

K 

k=l 

ch"'+^ 



\pk,l\ 



k<e. 



Pkeh-PkAlH-i^r":') - ^"^^ WWkWm+r^in") + WuiWrn+m^a")) 



c/ll+"(log/^)^(™)||pfe,,||^.+„^^,„ 
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where c is a constant independent of h and a, and (j{m) = if m < M — 2 and 
j3(m) = 1 ifm^M-l. 

Proof of lemma |4j We first give tlie proof for Pi finite element. The general 
proof is quite technical and is based on Lemma [H given in Appendix |B] Remind that 
we have denoted as x^'^ ,Xi'' , ...,xll'^i,xli^ the vertices of the triangulation of T^'^ 
that belong to F^'*^. To any r/^^^ in yf^'' n Hq{T^'^) we then associate the element 
in Wf^^ as follows 



ril,k\x, —Xr.' ) 1 l,k l,kr 

, t,k" over \xq' , x{ 

1 l.k e,k r 

r]i^k over \x^ ,a;„_i[ 
l^^pn^) over Jx„_i,<''=[ 



recalling that all norms are equivalent over the space of polynomials of degree 1 we 
easily deduce that there exists a constant c such that 

IIV'^,/c||L2(rk,f) < c||77£,fc||i2(rfc,f). 

Moreover, it is straightforward to derive 

{m,k + T^k,i{m,k))'4>i,k ^ I m.k'^i>,k+ / {TTk,e{m,k) f 



T^k.eive,k)iiJi,k - Vi.k) 



Then, by using the relation 



T^k,e{rii,k){'ip£,k - rjt.k) > -]^{TrkAVi,k)f - ^{i^e,k - m,kY 

leads to 

{ve,k + Trk,e{ve,k))'4'i,k > [ Ve,kiJe,k + 7; [ {T^k,i{m,k)f -\ ( {'4'i,k-m,kf 
We realize now that, over the first interval. 



Mi,k--{rP,,k-m,k))^ { ^,,k. --^ ,,,k ^^,kJ ^^.k 

J\x^ ,x{ [ - Xo j ^ \Xi - x„ ) 

noticing that 

[X — Xq ) 1 [X — X^ ) I [X — Xq ) 



1^0 ,^1 [ K^l — ^0 ) \Xl — Xq ) J]xg ,Xj_ I (Xi — Xq ) 



l,k i.k, , ,l,k fh. 



by recalling that ipi^k is constant on ]a;Q' ,Xi' [ and ]a:„'_i,2:„' [, we derive that 



ive.k + Trk,iive,k)He,k > Ve k 

Yk,e Jy''''^ 



which ends the proof of lemma 3) 
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PROOF of Lemma [4] in the general case Using the definition of TTfc ^, (|2.15p . 
it is straightforward to derive 

+ 'n-k,£{r]e,k){ipe,k - Vi,k)- 

Then, using the relation 

Trkji'ne,k){ipi,k ~ VLk) > -{T^kAvi^k))"^ - ^(V-'i?,fc - VLkf' 



leads to 

[m,k + 'n-kAvi,k))'4'Lk > [ m,kfpi.k - 7 / (V'^,fc - ve,k f- 

Remind that we have denoted as Xq'^jX^''^, ...,x^'^^,x'^'' the vertices of the triangu- 
lation of F^''"' that belong to F^^'^. By Lemma |5] of appendix B, and an easy scaling 
argument, there exists c, C > 0, ipi from [xq a;^''^] into M and V'n from [x^n-iT^n'^] 
into R such that 

/ / £,k\ / ^,/c\ / / £,k \ { i.k \ 1 

■)pi[x{ ) = r]{x{ ), -0Af(a;„'_i) = V(Xn^i) and 

1 1 r-i'" r-ii' 

Taking ^/j^^fe in W^^''' as follows 



Ve,k = < over \x{ 



, 1 e.k i,k 

ipi over Jxg , Xi 

l,k I, 

V-n over ]x„'_i,<''=[ 
proves LemmaHlwith ci = min(l, c) and C2 = max(l, C). 

Proof of theorem |4j In order to prove this theorem, let us build an element that 
will belong to the discrete space and will be as close as the expected error to the 
solution. Let be the unique element of Xj^ defined as follows : 

• i'^kh)\dn'' is the best fit of Uk over dil^ in 3^^' : 

• u\f^ at the inner nodes of the triangulation (in H,'') coincide with the interpo- 
late of Uk- 

Then, it satisfies 

\\ulh - Ukh^dn") < cft.^+"||ufc||^2+™(nfc), (3.8) 
from which we deduce that 

W^kh ~ ukh^n^^} + - UkWH^n") < ch'^^"'\\uk\\H'^+,^(^nk), (3.9) 
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and, from Aubin-Nitsche estimate 

114, ~ < c;i2+™hfclk2+„(o.). (3.10) 

We define then separately the best fit p^^, oipk^e = over each F*^'^ in Wj^'^. These 
elements satisfy for < m < Af — 1 the error estimate 

Wplth - PkAlLHr^^n ^ c/i^+"bfe,£||^i+„^j,,_,^ (3.11) 
IIpL/. -PMll^-i(p.,.) < ch'^n\PkA\Hh^^^r^,,y (3.12) 

But there is very few chance that {u^,p^J satisfy the coupling condition (|2.9p . This 
element of (^JlfcLi ^hj ^ (llt^i misses (|2.9p of elements ek,i and rji^k such that 

(Pkih + e/c,^ + aulMk.i = / {~p\kh + otm,k + au^/j)-^^/, ^^k,i S W^/J"'^ (3.13) 



(rife/^ + + aul,,)V'^,fc = / (-pLh - ^ki + auliJ^i^k, VV'f.fc e W';^'''. (3.14) 

In order to correct that, without polluting (|3.8[) - (|3.12[) . for each couple {k,£) we 
choose one side, say the smaller indexed one, hereafter we shall also assume that each 
couple {k,£) is ordered by fc < ^. Associated to that choice, we define ek^e € ^h'^i 
Vi,k e X''' n i?o^(r'='^), such that iuh,pj satisfy ([21]) where we define 



ueh = u\h + y^^Tle,k{r]i,k) 
k<e 



Pkih = plih + <^k,e (for k < £) (3.15) 

where TZg^k is a discrete lifting operator (see [36], [7j) that to any element of J^f'*^ fl 
7?g(r'^'^) associates a finite element function over that vanishes over dQ^\r'''^ and 
satisfies 

\\'RiAw)\\H^(n'-)<c\\w\\ ^ (3.16) 

where c is /i-independent. 

The set of equations (|3.13I1 - (I3.14|) for e^j and 77^^^ results in a square system of linear 
algebraic equations that can be written as follows 



{ek,i - Q:f]i^kHk,£ = / eiV'fc,^, VV'/c.f e W^'^ (3.17) 



{€k,£ + arii^k)iJe,k = / e2-0£,fc, VV'f^fc e Vy",^''' (3.18) 
pfc,^ jpfc,^ 
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with 

ei = -Plih - P\kh + "(wk - ^kh) (3-19) 

and 

62 = -plih - Pikh + a("fe/i - (3.20) 

Proposition 3.1. The linear system Ii3.17\ )- ^3J8\) is well posed. 
Proof: With these notations, (|3.17p yields 

ek,e ^ T^k,i{arie^k + ei) (3-21) 

and |37T8)) yields 

arji^k = T^t.ki-f-k^i + 62). (3.22) 

As (|3.17p - (l3.18p is a square linear system, it suffices to prove uniqueness for ei and 
62 null. From ((X^ - iP^ . we get 

= ?7£,fc + T^l,kT^kAVl,k)- 

SO that for all i/'£,fc in W^'^ , 



= / {m,k + '^k.t{m,k))^i,k- 

By Lemma 131 this proves that jy^^fe is zero, thus by (I3.2ip . is zero. 
Let us resume the proof of theorem 0) By (I3.2ip and (|3.22p we have 

[riLk + 'n-kj{'rie,k))ipeM = - / (e2 - TTk,e{ei))ipe,k, yipe,k S Wf^'''- (3.23) 

In order to estimate Hpfe^^ - Pfc^dl^-i^pfc.;-, and \\ueh ~ ue\\Hi(fi'), we first estimate 

\\''U,k\\L2(r''-i) ■ 

from (HH) and (|3:23l) we get 

ci||?7/!,fc||i2(rfc.«) < -||e2 - 7rfe/(ei)|ji2(rf=,<!)||V'£,fc||L2(r<=,f) (3.24) 



and using dSJ]) in ((3?24)) 



hence 



Ih^.fclUa/rfc.M < ^-||e2 - 7rfe.£(ei)||i2(r)=,f) 
^ ' aci ^ ' 



lhf,fc||L2(rfc,<!) < — (||e2||L2(rfc,*) + ||ei||i2(rfc,i)) (3.25) 



Now, from (|3l9l) and dS^Ql), for i = 1, 2 

||ei||i2(rfc,<!) < IbL/i +P£fc/JlL=^(r'=.f) + a||"k - '"LJlL2(r'=.0 
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and recalling that p^ i = — —-^ — —Pi,k over each V'^'^ 

WPkih +P£fc/JU2(rfc.o < \\pleh -Pfe/llL2(r'=.o + \\Pekh - PeM\L^{r'''') 
\Wih ~ '"willL2(r'=.0 < \\ulh ~ Ufe|lL2(r(=,«) + \\ulf^ - u£||^2(rfc,<!) 

so that, using (|3.8p and (|3.1ip . we derive for i — 1,2 and < m < M — 1 

||e.||i2(r.,.) < ca/ii+"(||ufe||ff2+„.(j,.) + ||u,||ff2+„(oo) + c/i^+"|bfc,f (3.26) 
and (f3?25l) yields for < m < M - 1 

3 chi~^"^ 



a 



We can now evaluate \\pkih - Pfe^^H^- i ^p^^^, using (|3.15|) 



WPkm -PkAH-iir".') - ll^'=^^llff-^(r'=.o + IbLh - Pfc.^ll^-ijr^,.)- (3-28) 

The term WpIhi — Pfe,f|| i , is estimated in p.l2|) . so let us focus on the term 
l|efe,dl^-i(pfc,i)- From (|3.2ip we have, 

||efc/||^_i^^^ ^^ < a||77^,fc||^_i^^^_,^ + l|ei||^-i^p, ,^ + \\{Id ~ Trk,i){aT]i,k + ei)\\ (3.29) 

To evaluate we proceed as for ||ei||i2(r'=.'!) and from p.lOp and p.l2p 

we have, for i = 1, 2 and < m < M: 

||e.||^-^(r.,,) < cah'+'^^iWukWH^+^^in^) + WwU^+^Hnn) + ^^^'^"'WPkA ^i-^^f^^.^ey (3-30) 
The third term in the right member of p.29p satisfies 

\\{Id- nkj){ave,k + ei)\\^_i^^^^^^ < cVh\\aTje^k + ei|li2(r)=,f). 
Then, using ([XTT]) and ((5:^ yields 

\\{Id- TrkA{c^Ve,k + eOll^-i^p^.^^ < cah'^+'^{\\uk\\ n^+^in") + IWeWn^+^in')) 

In order to estimate the term ||77^.fc|| „- 1 in (|3.29p . we use p.23p : 

H 2 (r '• ■) 

2 / m^k^i^k = 2 / T]i,kil'i,k 

-(/ {ve,k + TTk.e{ve,k))tpe.k - - I {e2 - T^k/ieiWi^k) 



that is 



2 / m,kipe.k ^ / {m,k - ■^k,e'rii,k)ipi,k -\ / (e2 - 7rfe.^(ei ))?/'<;, fc- 
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Using the symmetry of the operator tt^ ^ we deduce 

2 / r]e,ktpe,k = / {ipi.k - T!-k.eipi.k)r]£.k + - (e2 - 7rfe,£(ei))V'<!,fc- 
then, from ()2.16|) yields 

I J^^ ^ m,ki^e,k\ < cVh\\Tji^k\\mr''-')\\^e,k\\Hi(r>',i) + ^"^2 - TTfc.Kei)!!^- 1 (p,,,-,||V'Afcll^i(rk.^) 
and thus, we have 

ll%fc|lff-i(r.,.) ^ c\/)I|h£,fc||i2(r'».o + ^||e2 - 7rfc,Kei)||^_i ^^^^ 
Then, using (|3.27p and the fact that 

Ilea - nkAei)\\j^-i^^,,,^ < I|e2||^- i ^j,,,,, + ||ei|l^_ i + \\ei - nkj{ei)\\ ^^^^^^ 
with and ([3:301) yields for < m < M - 1 

Using the previous inequality in (I3.29p . p.28p yields 

\\Pkih-PkAH-i(^r>'.') - + ll^dk^+"(n^)) + c/i'+"l|Pfc.dl^^+™(r-.,)- (3-31) 

Let us now estimate \\ugh ~ ■ 

\\u£h ~ UiWH^n^) < - ue\\Hi{ir) + ^ \\T^£,k{ve,k)\\H^ni^) (3-32) 

fe<^ 

and from p.l6p 

\\T^Lk{m,k)\\m{n') < ch^.fell i 

then, with an inverse inequality 

\\T^i,k{rieA\\H^in>') < ch^^WmMlL^ir''.'!)- 
Hence, from (I3.27P we have for < m < M — 1 

\\T^eMm,k)\\min') < ch^^"'{\\uk\\H^+,,.^nk^ + \\ue\\H2+m^nn) 

+ — l|PMll^i+-,.(r.,.) 

and (|332)) yields 

\\u£h - U^|l_f/i(0«) < ch^^"'\\ui\\H2+n.^Qe^ + ch^+"' ^ ||ufe||^2+™(ofc) 

+— EllPMll^i+™(r.,.)- (3.33) 
fc<« 
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As M e H^+"\VL^), 



and for < m < M - 1 



\\Uh -U\\*< Y\ WuWm + r^in")- 

a ^ — ' 

z— 1 



Proof of theorem [5} The proof is the same that for theorem [H except that the 
relation (|3.1ip for < m < M — 1 is changed using the foUowing lemma 
Lemma 5. pj,^^ satisfy for Q < m < M — \ the error estimate 

\\pliH^PkA\L^-ir^-^) < ch^^"' (log/z)^(™)|bfc,,||^3^„^^, (3.34) 

Therefore, p.l2p is changed in 

and p.3ip is changed in 

WPkeh --PA;,dl^-i(pfc,f) < ca/i^+"(||ufc||H2+™(fi»=) + |Ih2+™(o<!)) 

W+"(log/i)'='(™)|b,,,||^3,„^^^,^ 

and (|3.33p is changed in 

\\ueh - ue\\Hi(<.v) < c/i^+"||u£||/f2+™(nf) + c/i^+" ^ ||wfc||^2+™(nfc) 

Proof of lemma [5} For < m < M ~ 1, j3{m) — and the estimate p.34p is 
standard. For to = M — 1, let f>kth be the unique element of W^'^ defined as follows : 

• (Vkih)\\ f.fc 1 coincide with the interpolate of degree M of pk i- 

• (Pfctt)|[3.^,fc .jf.fe] and {vkih)\y^>-,k ^ coincide with the interpolate of degree 
M - 1 of pk,i- 

Then, we have 

\\Vkih -Pfe,f|li2(rfc,*) < \\Vkih - Vk,i\\\'^{T^.ty 
Using Dcny-Lions theorem we have 

Let (y3 = '^-^r- In order to analyse the two extreme contributions, we use Holder's 
inequality : 

eh ^ 

(^2 < /l^"p||(/3||ip(o^;,). 
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Then, we use the estimate 

M 



where c is a constant. Thus we have 

Jo ' ~ ^ Hi{0,h) 



ip^ < cp'^h^^p\\(pf 1 



Then, 



£ k 



Now we take p — log h and thus we obtain 



'"'ij. 1^1' ^ '°e')^ii-'iii4"'rt-rf'.- 

In a same way we have 



,2M / i" ^fe,^ |2 / ,„„ j.\2||„_ j|2 



1^1^ s '"f')'"''""^*",,;^,.;-,' 



and thus we obtain 

\\VMh-VKl\\h-^{V^^t) < c(^*^+^log/l)||pfc,£|l^i+M(pfc,,y 

which ends the proof of lemma 

3.3. Error Estimates.. Thanks to p.Sp . we have the following error estimates: 
Theorem 6. Assume that the solution u of Il2.1\) - i2.^) is in i7^(ri) ni?g (r2), and 

Uk = U|f2fc e H'^+"'{n''), with Af - 1 > m > 0, and let pkj = over each V'"''^ . 
Then, there exists a constant c independent of h and a such that 



in") 



7 ' 

k=l I 



fc=l 

K 



Theorem 7. Assume that the solution u of h2.1]) - h2.A^} is in H^{n) n H^{Q), 
Uk = u^n" e H'^+"'{n''), and pkj ^ £^ is in H^+"'{rk,i) with M - 1 > m > 0. 
Then there exists a constant c independent of h and a such that 

K 

- Till , < r.Crv;)2+'" -J- 



fc=i 

K 



k=l i 
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with P{m) — if m < M — 2 and f3{m) = 1 if m — M — 1. 

Remark 1. Let us consider a Pi finite element approximation. If the solution 
u of i2.1\) - [K^) is in Yl^=i -ff* (^'^) o.iT-d a is a constant independent of h then, from 
theorem 

K 

\\uh-u\\* < WnWH^n") 

i=l 

and this result is not optimal. In order to improve this fact, we have to choose a 
parameter a which depends on h, or assume thatue Uk=iHU^''), 

or assume that 

ll^Uk^iHU^") andpk,, = 1^ € H'2{Tkj): 

• // the solution u of \2.1]) - \2.^) is in JlfeLi HI{VL^) and a = f , then 

K 

i=l 

• If the solution u of (Ml\l-(K^ is in Hf^i Hl{n'^), pk,i = is in iji {Tk,t) 
and a is a constant independent of h then 

\\uK-u\\* = 0{h\\og{h)\). 

3.4. Analysis of the best fit in 3D. In this section, we prove lemma |3] for a 
Pi-discretization in 3D. We shall use the construction proposed in [TT]. In order to 
make the reading easy, we shall recall the notations of the above mentioned paper. 
The analysis is done on one subdomain Vt^ that will be fixed in what follows. A typical 
interface between this subdomain and a generic subdomain 17/ will be denoted by T. 
We denote by T the restriction to T of the triangulation . Let S{T) denote the 
space of piecewise linear functions with respect to T which are continuous on F and 
vanish on its boundary. The space of the Lagrange multipliers on F, defined below, 
will be denoted by M(T). In 2D, the requirement dim M{T) =dim S{T) can be 
satisfied by lowering the degree of the finite elements on the intervals next to the end 
points of the interface. In 3D, it is slightly more complex (see [5]). Thus, we shall use 
the construction proposed in [11] in the case where all the vertices of the boundary 
of F are connected to zero or two vertices in the interior of F (figure [XT]) . Let V, Vo, 
dV denote respectively the set of all the vertices of T, the vertices in the interior of 
F, and the vertices on the boundary of F. The finite element basis functions will be 
denoted by <i>a, a e V. Thus, 

S{r) ^ span {$a : a e Vo}. 

For a e V, let CTa denote the support of $a, 

|J{^ e r : a e T}, 

and let Ma be the set of neighboring vertices in Vq of a: 

Wa ;= {6 e Vo : 6 e <7a}- 

Thus, 

u -^^ 

aedV 
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is the set of those interior vertices which have a neighbor on the boundary of F. If 
some triangle T e T has aU its vertices on the boundary of F, then there exists one 
(corner) vertex which has no neighbor in Vq. Let %, be the set of triangles T ^ T 
which have all their vertices on the boundary of F. For T € Tc, we denote by ct the 
only vertex of T that has no interior neighbour (such a vertex is unique as soon as 
the triangulation is fine enough) . Let Afc denote the vertices ar of M which belong 
to a triangle adjacent to a triangle T € 7^. Now, we define the space M{T) by 

M{T) := span aeVo}, 

where the basis functions are defined as follows : 



r $a, aeVo\Af 

bedVnaa 

+ Ab,aT'^b + ^CT a ^ are Afc 

fteavncTa^ 



the weights A^^a being defined in p.35p . For all boundary nodes c G dV connected 
to two interior nodes a and b, if Ta (resp. Tf,) denote the triangle having an edge on 
dV and a (resp. b) as the opposite vertex, then the weights are defined such that (see 

m) 

Ac.a + Ac,b = 1 and \Tb\Ac.a = \Ta\Ac,b- (3.35) 

M{T) is the notation introduced in [11], that we use here for the sake of clarity. Cor- 
responding to our previous notation, M{T) = W^'^ . 

To any u G S{T), u — X^aeVo we associates G M{T) where w = X^aeVn '"(^)^a- 

More explicitly, that means that to any u E S{T), we associate an element v E M{T) 
as follows (see figure [3TT|) : 

(i) V \s a. piecewise linear finite element on T 

(ii) for all interior nodes a, v{a) :— u{a) 

(iii) for all boundary nodes c, by assumption we have two situations: 

- c is connected to two interior nodes denoted by a and b. 
Then, v{c) :— Au{a) + Bu{h) where 

A + B =1 imA\Tb\A=\Ta\B (3.36) 

where Ta (resp. Tb) is the triangle having an edge on dV and a (resp. b) as the 
opposite vertex. 

- c is not connected to any interior point. We consider the triangle adjacent 
to the triangle to which c belongs to. This triangle has one interior node denoted 
by b. Then, we define v{h) :— u{h). 

Then, using the uniform regularity of T, it is easy to check that there exists a constant 
c independent of h such that 

M\l^t) < c||u|lL2(r). 



We shall need the following technical assumption: 

Assumption Let < C < 2/3. For any triangle Tc' having all three vertices on the 
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boundary ofT (see fiaure \3 . 1]] . we consider the two triangles ^ andTj/ ,, surrounding 
Tc' ■ We assume that 



1 C 

— min(|r,,c|, \Tj',c\) > 



In order to prove lemma [31 we prove the following lemma: 

Lemma 6. Let < C < 2/3, we assume the above assumption and that T is 
uniformly regular. Let u G S{T) and let v € M{T) constructed from u as explained 
above ( (i)-(iii)). 

Then, there exists c > such that, 

j {uv~^{u-vf)>c j u^. (3.37) 

Proof of lemma [6j Let us introduce the notation 

C 

{uv - —{u - vf). 



We have 

QT = \j{u + vf - (1 + 2C){u - vf 
In order to estimate Qr, we remark that 

TeT 

where 



Qt = \ j {u + vf-{\ + 2C){u - vf. 
It 
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We have four kinds of triangles: 

1. Inner triangles i.e they don't touch the boundary of T . 

2. Triangles which have only one vertex on the boundary 

3. Triangles which have two vertices on the boundary 

4. Triangles which have three vertices on the boundary 

Inner triangles. On an inner triangle T, u = v so that for all C > 0, we have 

Qt>c 

Jt 

for c < 1. 

Triangles having only one vertex on the boundary. Let Tj ^ be such a triangle (see 
figure [5?T|) . We use the following notations: Ui = u{a) = v{a), Ui+i = u{b) = vib) and 
Vi = v{c). First notice that we have (remember u{c) — 0) 



= ^-^^ (u"^ + uf+i + (ui + Ui+if^ = -^^^ ( 2u- + 2uf^i + 2uiUi+i 
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see for example pTO] (II. 8. 4). As for Qy. ^, we have 



+ (2uj + 2uj+i + Au^ + Bu,+iY - 2(1 + 2C){Au^ + Bu,+i) 
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8uf + 8u^_^i + 8uiUi+i + 4(uj + Mi+i)(Au.i + Bm+i) 

-4C{Au, + Bu,+if 



If we take C = 1 and use ^ + _B = 1, we get: 
48 



Qt,,, = [ 4u- + 4m^+i + Au^U^+l + AAB{ui - U^+lY + 4(m, + U^+lf 



1 



Hence, for all < C < 1, we have: 



Therefore, 



for < C < 1 and < c < 1/2. Wc shall also use in the sequel the estimate: 

Qt.,. > (3.38) 

Triangles having two vertices on the boundary. We consider now a triangle Ti^r 
having two vertices on the boundary of the face F, see figure [3?11 Let Mr = {«, Ti^r 
has two vertices on the boundary of F}. First notice that we have 

2 ^ ^ ^2 



And we have 
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4u^_^i + w,^ + Wj^+i + (2ui+i + + Ui+i)^ 



(l + 2C)(i;?+i;2 +(„^+„,,+,)2) 



48 

Then, 



Qt„. > %^ ( 8^.2^1 - &Cv} - 6Cvf+, + 4u,+,{tH + ) . 
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Defining Ei := Wi+iWi and Fi := u^+iUj+i (cf. [TT] page 11), we have: 



Now we sum these terms over ah the triangles having two vertices on the boundary 

r. 

> / + ^ E (|r.,r|(-6CT;2 + + \T,^,,.\i-6Cv^ + m-i)). 

Ju.eM^T.,- 48 ^^^^ 

(3.39) 

The condition p.36p leads to the inequality 

\T,^r\E, + |T,_i,,|F,_i = {\T,^r\ + \T,^l,r\)vf 

(see equation after (3.19) in [11), so that we get: 

\T,^r\{-6Cvf + AEi) + |r,_i,,|(-6Cf2 + 4^,„i) = (|r,,,| + |r,_i,,|)(4 - &C)vl 
This term cancels for C = 2/3. Hence for < C < 2/3, inequality p.39|) becomes: 



E Qt... > / 



Therefore, for < C < 2/3 and < c < 1, 

E Qt^. ^ c / 



Triangles having all three vertices on the boundary. Let Tc' be such a triangle 
(see figure IXT|) . We have to control: 

by the integrals over the two triangles Ti^c and Tji^c surrounding Tc'. This can be 
achieved using the assumption 

-^min(|r,,,|,|T,.,,|)>^|r,n 

and using that fi'om (|3.38[) . we have 



In conclusion, we have that (|3.37p holds with c = 1/4 for < C < 2/3. 
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^ Boundary of T 



Fig. 3.1. Two different situations of 2D triangulation of the interface F, next to it's boundary 
(near cross points) 

4. Numerical results. On the unit square ft = (0, 1) x (0, 1) we consider the 
problem 

{Id - A)w(a;, y) ^ ^^(y^ - 2) - Gxy^ + (i + _^ y^)sin{xy), {x, y) e ^, 
u — x^y^ + sin{xy), {x, y) G 90, 

whose exact solution is u(x,y) = x^y^ + sin{xy). We decompose the unit square 
into non-overlapping subdomains with meshes generated in an independent manner. 
The computed solution is the solution at convergence of the discrete algorithm (|2.1ip - 
(|2.12p . with a stopping criterion on the jumps of interface conditions that must be 
smaller than 10^^. 

Remark 2. In the implementation of the method, the main difficulty lies in 
computing projections between non matching grids. In US}/ we present an efficient 
algorithm in two dimensions to perform the required projections between arbitrary 
grids, in the same spirit as in '191 for finite volume discretization with projections on 
piecewise constant functions. 

4.1. Choice of the Robin parameter a. In our simulations the Robin pa- 
rameter is either an arbitrary constant or is obtained by minimizing the convergence 
rate (and depend of the mesh size in that case). In the conforming two subdomains 
case, with constant mesh size h, the optimal theoretical value of a which minimizes 
the convergence rate at the continuous level is : 

In the non-conforming case, has the mesh size is different for each side of the interface, 
we consider the following values : 

a™„- [(7r2 + l)((-^)2 + l)]3 

'^mean 

where hmin, hmean and hmax stands respectively for the smallest meanest or highest 
step size on the interface. 
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4.2. An example of computed solution. We decompose the unit square into 
four non-overlapping subdomains with meshes generated as shown in Figure B?T1 The 
Robin parameter is a = 10. On Figure 14.21 we show that the computed solution is 
close to the continuous solution. 



Mesh 




D 0.2 0.4 0.6 0.8 1 



Fig. 4.1. Domain decomposition with non- conforming meshes. 



Computed solution 




Fig. 4.2. Computed solution. 

4.3. error between the continuous and discrete solutions. In this 
part, we compare the relative error in the non-conforming case to the error ob- 
tained on a uniform conforming grid. 

Definition of the relative error : Let K be the number of subdomains. Let 
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Ui = M|f2i, I < i < K (where u is the continuous solution), and let = (u/i)|n» 

where Uf^ is the solution of the discrete problem (|2.10p . Now, let E^r^ — ||w||, and let 
El = \\{uh)i - ■"illHi(O'), l<i< K. Let 

i=l 

The relative error is then E/Eex- 

In this example, we take for the Robin parameter. This choice is motivated 

by the results of section but we obtain similar results in the case of a = 10. 
We consider four initial meshes : the two uniform conforming meshes (mesh 1 and 
4) of figure 14.31 and the two non-conforming meshes (mesh 2 and 3) of figure 14.41 
In the non-conforming case, the unit square is decomposed into four non-overlapping 
subdomains numbered as in figure l4?5l 




Fig. 4.3. Uniform conforming meshes : mesh 1 (on the left), and mesh 4 (on the right) 




Fig. 4.5. N on- overlapping domain decomposition of the unit square 
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Figure [321 shows the relative error versus the number of refinement for these 
four meshes, and the mesh size h versus the number of refinement, in logarithmic 
scale. At each refinement, the mesh size is divided by two. The results of figure 14.61 
show that the relative error tends to zero at the same rate than the mesh size, 
and this fits with the theoretical error estimates of theorem |4l On the other hand, 
we observe that the two curves corresponding to the non-conforming meshes (mesh 2 
and mesh 3) are between the curves of the conforming meshes (mesh 1 and mesh 4). 
The relative error for mesh 2 is smaller than the one corresponding to mesh 3, 
and this is because mesh 2 is more refined than mesh 3 in subdomain f2*, where the 
solution steeply varies. 



HI relative error versus the number of refinement 




number ot refinements 



Fig. 4.6. relative error versus the number of refinements for the initial meshes : mesh 1, 
(diamond line), mesh 2 (solid line), mesh 3 (dashed line), and mesh ^ (star line). The triangle line 
is the mesh size h versus the number of refinements, in logarithmic scale 



More precisely, let us compare for mesh 2, the relative error in the domain 
Vt^\jVL^\J Vt^ to the relative H'^ error in the subdomain Q!^ (which is the subdomain 
where the solution steeply varies). This comparison can be done in Table 1. 



Refinement 


{El + El + Elf/^/E,^ 


E4,/Eex 


E/Eex 





1.45e-01 


1.46e-01 


2.06e-01 


1 


7.17e-02 


7.02e-02 


1.004C-01 


2 


3.59e-02 


3.49e-02 


5.01e-02 


3 


1.79e-02 


1.73e-02 


2.49e-02 


4 


8.73e-03 


8.46e-03 


1.21e-02 



Table 1: Comparison, in the case of mesh 2, for different refinements (column one), 
of the relative error in the domain composed by subdomains , Vl^ and 51'^ 
(column 2) to the relative error in the subdomain fi'* (column 3). The fourth 
column is the relative error in the whole domain. 
We observe that, as expected, the relative error in the domain composed by 
subdomains fi^, fi^ and (second column of table 1) is close to the relative 
error in the subdomain fi^ (third column of table 1). Indeed, the mesh 2 is more 
refined in the subdomain il^ where the solution steeply varies. 
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Refinement 




Ei/Eex 


E/Eex 





1.26e-01 


2.04e-01 


2.40C-01 


1 


5.57e-02 


1.04e-01 


1.18e-01 


2 


2.74e-02 


5.22e-02 


5.90e-02 


3 


1.36e-02 


2.59e-02 


2.93e-02 


4 


6.64e-03 


1.26e-02 


1.43e-02 



Table 2: Comparison, in the case of mesh 3, for different refinements (column one), 
of the relative error in the domain composed by subdomains Sl^, Vl^ and Vl^ 
(column 2) to the relative error in the subdomain Vl'^ (column 3). The fourth 
column is the relative error in the whole domain. 



Let us now do the same comparison in the case of mesh 3. This mesh is coarser in 
the subdomain fi'* where the solution steeply varies. In table 2, we observe that as 
expected, the relative error in the domain composed by subdomains f2^, Vl^ and 
fi'^ (second column of table 2) is smaller (almost half) than the relative error in 
the subdomain fi"* (third column of table 2). That one is close to the relative 
error in the whole domain (fourth column of table 2), because mesh 3 is coarser in 
the subdomain where the solution steeply varies. 



4.4. Convergence : Choice of the Robin parameter. Let us now study 
the convergence speed to reach the discrete solution, for different values of the Robin 
parameter a. We first consider a domain decomposition in two subdomains, and then 
in four subdomains. 



AAA. 2 subdomain case. In this part, the unit square is decomposed in two 
subdomains with non-conforming meshes (with 81 and 153 nodes respectively) as 
shown in figure 14.71 On figure 14.81 we represent the relative error between the 
discrete Schwarz converged solution and the iterate solution, for different values of 
the Robin parameter a. We observe that the optimal numerical value of the Robin 
parameter is close to amean and near and amax- 
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Mesh 




Fig. 4.7. Domain decomposition in 2 subdomains with non-conforming grids 



As the relative H error didn't show where the error is highest, we also look at the 
relative L°° error between the discrete Schwarz converged solution and the solution 
at iteration p, for different values of the Robin parameter a. We obtain similar results 
as for the relative error (see figure l4!9|) . 

The Schwarz algorithm can be interpreted as a Jacobi algorithm applied to an interface 
problem (see [32]). In order to accelerate the convergence, we can replace the Jacobi 
algorithm by a Gmres CJS]) algorithm. Figures [4.101 and 14.111 show respectively the 
relative error and the relative error between the discrete Gmres converged 
solution and the iterate solution, for different values of the Robin parameter a. In the 
case where a = ameam we observe that the convergence is accelerated by a factor 2 for 
Gmres, compared to Schwarz algorithm. Also, the gap between the error values for 
different a is decreasing when using Gmres algorithm, compared to Schwarz method. 
The Gmres algorithm is less sensitive to the choice of the Robin parameter. 
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Schwarz ; relative error versus the number of iterations, for different vaiues of a 




20 40 60 80 100 120 140 160 

number of iterations 



Fig. 4.8. relative error between the discrete Schwarz converged solution and the iterate 
solution, for different values of the Robin parameter a 



Schwarz : reiative C error versus the number of iterations, for different values of a 




20 40 60 80 100 120 140 160 

number of iterations 



Fig. 4.9. Relative L°° error between the discrete Schwarz converged solution and the iterate 
solution, for different values of the Robin parameter a 
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Gmres ; relative error versus the number of iterations, for different values of a 




2 4 6 8 10 12 14 16 18 20 

number of iterations 

Fig. 4.10. relative error between the discrete Gmres converged solution and the iterate 
solution, for different values of the Robin parameter a 



Gmres : relative L°° error versus the number of iterations, for different values of a 




2 4 6 8 10 12 14 16 18 

number of iterations 



Fig. 4.11. Relative L°° error between the discrete Gmres converged solution and the iterate 
solution, for different values of the Robin parameter a 

4.4.2. 4 subdomain case. In this part, the unit square is decomposed into four 
subdomains with non-conforming meshes (with 189, 81, 45 and 153 nodes respectively) 
as shown in figure HTT^ On figure and HTT^ respectively, we represent the relative 

error and the relative L°° error between the discrete Schwarz converged solution 
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and the iterate solution, for different values of the Robin parameter a. We observe 
that the optimal numerical value of the Robin parameter is close to amean a-nd near 
OLrain and cimax^ as in the two subdomain case. 




4.4.3. Conclusions. The numerical results on the relative error between 
the continuous and discrete solutions correspond to the theoretical error estimates of 
theorem|4l On the other hand, we observe that, for a fixed number of mesh points, the 
relative error between the continuous and discrete solutions is smaller for a mesh 
refined in the region of the domain where the solution steeply varies, than for a mesh 
which is coarser in that region. In term of convergence speed to reach the discrete 
solution, the Robin parameter a must depend of the mesh size, and our simulations 
show that a = amean is close to the optimal numerical value. 
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Schwarz : relative H error versus the number of iterations, for different vaiues of a 




Fig. 4.13. relative error between the discrete Schwarz converged solution and the iterate 
solution, for different values of the Robin parameter a 



Scfiwarz : relative error versus tfie number of iterations, for different vaiues of a 



10.6708 (a ) 
15.3466 (a ) 
28.7856 la^ . ) 
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60 80 1 00 

number of iterations 



120 140 



Fig. 4.14. Relative L°° error between the discrete Schwarz converged solution and the iterate 
solution, for different values of the Robin parameter a 



Appendix A. Inf-sup condition.. The purpose of this annex is to show that 
the proof of [Ij can be extended to the 3D situation. Indeed the main ingredients 
required for the extensions have been proven in pjj. Let us first recall a standard 
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stability result in higher norms of the projection operator ttj. £ from L'^(r'^'^) onto 



y^'^ n i?o(r'''0 orthogonal to W^'^ 



Lemma 7. Making the hypothesis that the triangulation is uniformly regular, 
there exists a constant c > such that 

Proof From (|3.37p we deduce a uniform inf-sup condition between y^'^ n Hq{T'''^) 
and Wl^'^ in L'^(r^'^). It results that the projection operator tt is stable in L^(r'^'^) 
and thus there exists a constant ci > such that 

V^; e H|o(r'''), 11^ - TTkMLHr^^n < ci/i^ Ikll i , • 

Let TTk^e denote the orthogonal projection operator from i?QQ(r'^'^) onto 3^,j'^ni?Q (F*^'^) 

1 1 
for iJQ^Q(F'^'^) inner product. Then, for all v in H^q{T'''^), 

\\Ttk fv\\ 1 <l|7rfc£w|| 1 + \\T^k eV — TTk iv\\ 1 

Then, with an inverse inequality, there exists a constant C2 > such that 

II^M^II^I^j,,,,^ - Il"ll4(r...) +'^2/^-^||^M^^-^M«||L^(^'=.*)• 



Thus, 



iTrfc^wll 1 < llwll 1 +C2/1 -^c'h^WvW 1 , 



and then, with c — 1 + c'c2 ,we have 

ll^M^'ll i <c|li'|| 1 ,Vvei/|o(r'=''), 

which ends the proof of lemma [T] 

— - ~ k i 

Then from the definition of the iJ* ^j-pfe.f-j norm, for any ph,k,i in Wf^' , there ex- 
ists an element w'^'^ in Hq^{T'''^) such that 



and w'''^ can be chosen such that 

We apply now the projection operator on w'''^ from lemma [71 We derive that 
TTkA^"'') = G 3^^' n Fo'(r'='0 and 

II W^^ll i < C||P/!. fc ^11 1 , 
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and 

Ph,k,ewy = / Ph.k,iw^'^ = WPh.k.A? 1 

It remains to lift w^^ over fi'^, this is done by prolongating lu^^ by zero over dVl^\T^'^ 
and lifting this element of H^{d^l'^) over il*^ as proposed in [7]. 

Appendix B. Extension in 2D for high order approximations. 

Lemma 8. Let 1 < p < be an integer. There exists c and C > such that for 
all r] e P*'([-l, 1]) s.t. r/(-l) = there exists ijj G FP'^{[~1, 1]) s.t. 



'?(!)= ^(1) 



and 



ilr <C r)\ 
J -I 

This lemma has been proven in the case p = 1 at section [321 For p > 2, we prove 
this lemma by studying for a given rj G P''([— 1, 1]), 77 the maximization problem 



Find V e PP" !([-!, 1]) such that 

J{'ip;ri) — max J{(f>;ri). (B.l) 

0ePf-i(hi,l]) 
0(1) = 77(1) 

The function J is strictly concave and there exists a function satisfying the constraint. 
This problem admits a solution. The functional J{4>, rj) being quadratic in (0, 77) and 
the constraint being afhne, the optimality condition shows that the problem reduces 
to a linear problem whose right hand side depends linearly of 77. The affine constraint 
being of rank one, the problem (jB.ip admits a unique solution which depends linearly 
of 77. Therefore, it makes sense to introduce the operator: 



77 i-^ -0 solution to (|B.l 



where Pg([— 1, 1]) is the set of functions of P''([— 1, 1]) that vanish at —1. In Lemma[SJ 
we take ip — S{ri). The operator S is linear from a finite dimensional space to another 
so that it is continuous for any norm on these spaces. Therefore there exists C > 
such that J^^ tA^ < C J^^ 77^. Moreover, the function 

i?:Pg([-l,l])\{0} ^ M 

77 
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is continuous and such that H{ri) = H{ari) for any a ^ 0. Therefore, it reaches its 
minimum and proving Lemma [5] amounts to prove 

Lemma 9. Let p < 13 and rj e PP([— 1, 1]) s.t. ?/(—!) = and 77 is not the null 
function. 
Then, 

J{S{rj);7j)>0. 



Proof. We make use of the Legendre polynomials 

Lo{x) = 1, Li{x) = X, {m + l)L,n+iix) = {2m+ I) X Lm{x) ~ mLm-iix), m>l 
Let us recall that for any m > 0, 

L,n{l) = 1, im(-l) = (-1)^ 



1 2 

/_ ^ Ljn {x) Lm' {x) dx = Sr, 



Jrti rn' 



2m +1 

The polynomial rj is decomposed on the Legendre polynomials 

p 

and ip — S{ri) is sought in the form 

p-i 

m=0 

SO that it maximizes the quantity J{ip]ri) under the constraint 77(1) = ipi^)- This 
corresponds to the min-max problem 

max min u) 

^gpp-i([_i,i]) ^eB 

where 

Ci^,^i) = J{^P;r1)-^i{^i^)-r^{^)). 

We have to prove that the optimal value is positive. The optimality relations w.r.t ip 
give 

3, , 1 , 2m + 1 

-^(Vni + 1lm+lj - -^Vm = M ^ , 1 < m < p - 1 



and 



Therefore, we get 



3 1 , M 
2^^^ 2^° = 2 • 



tp = Si] - 3rip Lp - fiRp^i (B.2) 
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p-l 

where ~ (2m + l)Lm and Il^^p-i|li2(]_i = 2p^. Hence, the dual problem 

m=0 

writes 

min G(u: 77) 

where 

G(/i;7?) := J{3T]-3T]pLp - fiRp-i;r]) - ^(V'(l) - 
and "0 satisfies (jB.2p . After some calculations, we get 

G(m; v) = yM' - M(2r/(1) - 3ryp) + {2MI,q_,,,^^ - l^^)- (B.3) 

The leading coefficient of G{fi; rj) is positive so that proving min^ G{fi; 77) is positive 
(and hence Lemma 1^1) is equivalent to prove 

Lemma 10. For p < 13, the discriminant of liB.3\) : 

A(r;):=(2,y(l)-3r;p)2+p2(_4||,,||2.^^j_^^^j^+9_^) (5.4) 

is negative if rj € PP([— 1, 1]), 77(— 1) = and rj is not the null function. Proof. We 
first treat separately the case p ~ 2. In this direct computation shows that 

80 2 40 133 2 

The discriminant of the corresponding bilinear form is —8632/9. It is negative and 
the lemma is proved in this case. 

We consider now the case p > 3. Let us introduce the vector space = {?? <E 
PP([— 1,1]) s.t. 77(— 1) ~ 0}. The function A(77) is quadratic so that it suffices to 
study the extrema of A{ri)/\\ri\\'j^2Q_i iq over Qp or equivalently to prove that the 
associated symmetric quadratic form in negative, i.e. its eigenvalues are negative. 
They correspond to the Lagrange multiplier solutions /ii of the following min-max 
problem 

min max £e(?7, ui) (B.5) 

where 

^e{v,fJ-i) A(?7) - Mi(hlli2Q„i^i[) - 1). 
We have to prove that fii < 0. We have 

= 2(277(1) - 377p)(2<577(l) - 3,577^) + p^{-8 <fj,Sr,> +18^^^) - 2^1 < 77, Si] > 

where < , > denotes the scalar product on i^(] — 1, 1[) and Sr] e Q^. 

Let us consider the vector space (1 — x'^)PP^^ C Q^. Any function 7 in (1 — x^)PP~^ 
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satisfies 7(— 1) — 7(1) — and 7p — 0. The optimality relation w.r.t. to (1 — x^)V^ ^ 
gives 

(-8p2 - 2^1) < &r) 0, V(5?7 e (1 - a;2)pp-3. 

We have either /ii — — < or 77 solution to (|B.5[) belongs to the space {(1 — 
^2-jpp-3|± p pp^ rp-j^g gj.g^ gg^gg corresponds to a negative value for /ii which is in 
agreement with the lemma to be proved. Let us study the latter case. We shall make 
use of 

Lemma 11. 



y ^'m L'm' (1 - x^) dx = 0, m', 



(B.6) 



L'J =m{m + l), (B.7) 
1 

1 

-^m Kn+1 = 0, (B.8) 

-1 



m(TO-l), (B.9) 



see II]/. From Lemma [TTl it can be proved that 
Lemma 12. 

{(1 - a:2)pf-3}^ n PP = Span{Lp, L'^, L'^_^}. 
Proof. From (jB.6p . it can be checked easily that 

{(1 - a;2)pp-3}± n PP = Span{i;+i, L;, L;_J. 
Moreover, we have 

= (2p + i)ip(x) + i;_i(x)) 

and thus lemma [T^ 

Therefore, there exists Ai,A2, A3 € M s.t. -q — XiLp + X2L'p + AaLp.;^. Since 77 is 
defined up to a constant and we only have to consider the two cases Ai = 1 or Ai = 0. 
Case 1 Ai = 1 
From 1) = 0, we get 

1_A,^^ + A3fc^=0, 



so that 



A2 = ^— + A3^. 

p{p +1) p + 1 



^ ^ {p-l)p^ (P^ + 1) ^^2 _ (24p4,20p3_8p^ + 4p) 
p+1 ^ (p+l)(2p+l) 

29 p2 + 13p- 1 - 



(p+l)(2p+l) 
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Since p is supposed larger than 1, the leading coefficient of A (77) is negative. If 
the discriminant of A(t/) is negative, the polynomial is negative for any A2. This 
discriminant has the value 



^g(pl_-13p-8)(p-l)p3 



2p+l 



and is negative for 2 < p < 13. 

Case 2 Ai = 
Prom 1) = 0, we get 



so that 



■a /(p + i) +A3 ^(p-i) =o, 



A2 = As^ 



p+1 

Since 77 is an eigenvalue, it is not zero and the above relation shows that we can take 
A3 = 1. Then, we have A2 = so that 
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Computed solution 




Gmres : relative H error versus the number of iterations, for different values of a 
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Gmres : relative L°° error versus the number of iterations, for different values of a 
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